t 


FI^TR-93-2166 

EmrlrooMental  Research  Papers,  No.  1126 


AD-A276  690 

lillllllll 


SENsmvmr  analysis  of  the  avhrr 

INFRARED  CIRRUS  MODEL 


Allan  J.  Buss^ 

Robert  P.  d'Entremont 


Approved  for  Public  Release;  distribution  unlimited 


PHILLIPS  LABORATORY 

Directorate  of  Geophysics 

AIR  FORCE  MATERIEL  COMMAND 

HANSCOM  AIR  FORCE  BASE,  MA  01731-3010 


*^4  *>  0*^ 


94-04232 

lililiilii 


Best 

Available 

Copy 


This  technicai  report  has  been  reviewed  and  is  approved  for  publication.' 


JC^EPM  W.SNOW 

Cmef.  Satellite  Meteorology  Branch 

Atmospheric  Sciences  Division 


ROBERT  A. /McCLATCHEY 
I^rector.  Atmospheric  Sciences  Division 


This  report  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technicai 
information  Center.  Ail  others  should  apply  to  the  National  Technicai  Information 
Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization,  please  notify  PU/TSI,  29 
Randolph  Street,  Hanscom  AFB,  MA  01731-3010.  This  will  assist  us  in  maintaining  a 
current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  requires  that  it  be  returned. 

'This  technical  report  has  been  reviewed  smd  is  approved  for  publication.' 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  coKection  of  information  is  estimated  to  atfera9e  i  hour  per  response,  including  the  time  for  reviewir>9  instructions,  searching  emstm^  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  coiieaion  of  information  Sertd  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
coltection  of  information,  including  suggestions  for  r^ucing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  information  Operations  and  Reports.  121$  Jefferson 
Davis  Highway.  Suite  1204.  Arlington,  VA  22202'4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704^)188).  Washington.  DC  20S03. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

19  July  1993  Scientific,  Interim 

4.  TITLE  AND  SUBTITLE 

Sensitivity  Analysis  of  the  AVHRR  Infrared 

Cirinis  Model 

5.  FUNDING  NUMBERS 

PE  62101F 

PR  6670 

TA  17 

WU  12 

6.  AUTHOR(S) 

Allan  J.  Bussey  and  Robert  P.  d'Entremont 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Phillips  Laboratoiry  (6PAS) 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

PLr-TR— 93-2166 

ERP.  Ho.  1126 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS{ES) 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  Public  Release 

Distribution  Unlimited 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

We  have  studied  the  sensitivity  to  errors  of  a  cirrus  cloud 
analysis  technique  using  NOAA  Advanced  Vexry  High  Resolution 

Radiometer  (AVHRR)  multispectral  imagery.  The  sensitivity 
analysis  reveals  possible  errors  in  confuting  optical  depths  and 
altitudes  for  nighttime  cirinis  cloud  data  samples  taken  from  the 

First  ISCCP  Regional  Experiment  (FIRE)  Phase  I  intensive  field 
obseirvations .  The  Appendix  contains  plots  of  the  computed 
errors  found  in  the  study. 

14.  SUBJECT  TERMS 

AVHRR  satellite  imagery.  Cirrus  cloud.  Cirrus 
altitude,  Emissivity,  Optical  depth 

15.  NUMBER  OF  PAGES 

42 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UHCLASSIFIED 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UNLIMITED 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std  Z39-18 


298-102 


Contents 


1.  INTRODUCTION  .  1 

2.  NOTATION .  2 

3.  CIRRUS  PROPERTY  ANALYSIS  METHOD  .  3 

4.  SENSITIVITY  ANALYSIS  METHOD .  12 

5.  SENSITIVITY  ANALYSIS .  13 

6.  SENSITIVITY  RESULTS  FOR  OPTICALLY  THIN  CIRRUS  CLOUDS: 

AVHRR  CHANNEL  4  EMISSIVITY  OF  0.2 .  15 

6.1  Surface  Temperature  of  277.1  K .  15 

6.2  Surface  Temperatvu*e  of  274.1  K .  16 

6.3  Surface  Temperature  of  271.1  K .  17 

7.  SENSITIVITY  RESULTS  FOR  OPTICALY  THIN  CIRRUS  CLOUDS: 

AVHRR  CHANNEL  4  EMISSIVITY  OF  0.5 .  18 

7.1  Surface  Temperature  277.1  K .  18 

7.2  Surface  Temperature  274.1  K .  18 

7.3  Surface  Temperature  271.1  K .  19 

8.  SENSITIVITY  RESULTS  FOR  OPTICALLY  THIN  CIRRUS 

CLOUDS:  AVHRR  CHANNEL  4  EMISSIVITY  OF  0.9 .  20 

8.1  Surface  Temperature  277.1  K  .  20 

8.2  Surface  Temperature  274.1  K .  21 

8.3  Surface  Temperature  271.1  K .  22 

9.  SUMMARY .  22 

BIBLIOGRAPHY .  25 

APPENDIX  A  .  26 


Aooession  For _ j  | 

■TIS  8PvA&I  ar 

DTIC  TAB  □ 

Just  If  lo"  ..loii _ 

1 

- - - 1 

Ry - i 

_ j 

Ava liability  .iOtse 
jAvaii  ai'd/cr 


iii 


Illustrations 


1.  AVHRR  10.7;un  Channel  4  Theoretical  Brightness  11 

Temperatures  as  a  Function  of  the  Channel  4  Emissivity 

^tie  Effective  Cirrus  Altitude 

2.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Hei^t  Errors  (km)  as  a 
Function  of  Brighi^ness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^  =  277.1  K  (AT^=  +3  K) 

and  True  Cirrus  Emissivity  84  =  0.2.  28 

3.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Hei^t  Errors  (km)  as  a 
Function  of  Briefness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^  =  274.1  K  (AT^=  0  K) 

and  True  Cirrus  Emissivity  84  =  0.2.  29 

4.  Plots  of  AIRC  •  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Height  Errors  Osm)  as  a 
Function  of  Briefness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^  =  271.1  K  (AT^=  -3  K) 

and  True  Cirrus  Emissivity  84  »  0.2.  30 

5.  Plots  of  AIRC  •  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^  =  277.1  K  (AT,fe=  +3  K) 

and  True  Cirrus  Emissivity  84  =  0.5.  31 

6.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Hei^t  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^  =  274.1  K  (AT,fc  =  0  K) 

and  True  Cirrus  Emissivity  84  =  0.5.  32 

7.  Plots  of  AIRC  •  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Hei^t  Errors  (km)  as  a 
Function  of  Briefness  Temperature  Measurement  Errors 
AT3  and  AT4.  Error  Plots  are  for  T^j.  =  271.1  K  (AT,fc=  -3  K) 

and  True  Cimis  Emissivity  84  =  0.5.  33 


iv 


8.  Plots  of  AIRC  -  Retrieved  Cirrus  Emisshrity  Errws 
(dimensionless)  and  Cloud  Hei^^t  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors 
AT,  and  AT4.  Error  Plots  are  for  =  277.1  K  (AT^  =  +3  K) 

and  True  Cirrus  Emissivity  =  0.9.  34 

9.  Plots  of  AIRC  •  Retrieved  Cirrus  Enussivity  Errors 
(dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Bri^^tness  Temperature  Measurement  Errors 
AT,  and  AT4.  Errca:  Plots  are  for  T,,^  s  274.1  K  (AT^  =  0  K) 

and  True  Cirrus  Emissivity  =  0.9.  35 

10.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors 
(dimensionless)  and  Cloud  Hei|^t  Ekitars  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors 
AT,  and  AT4.  Error  Plots  are  for  T^  =  271.1  K  (AT,fc=  -3  K) 

and  True  Cirrus  Emissivity  £4  =  0.9.  36 


V 


Tables 


1.  Corresponding  Values  of  Emissiviiies  for  5 

AVHRR  Channels  3  and  4,  listed  as  a  Function 

of  Cloud  Thickness. 

2.  Channels  3, 4,  and  5  Brightness  Temperature  8 

Measurements  for  Various  Effective  Cloud  Top 

Heists  and  Channel  4  Emissivities,  Calculated- 
Using  E!qs.  (1)  and  (4)  for  the  Madison, 

V^sconsin  Sample  Atmosphere. 

3.  £4  =  0.20;  CH  3  and  CH  4  AT  =  ±5  K  Where:  11 

T^  is  the  Surface  Temperature  (IQ;  is 

the  Effective  Cloud  Altitude  (km);  A  is 
the  Difference  (+  or  -)  in  Channel  4 
Emissivity  From  Expected  Values;  and  A  Z 
is  the  Difference  (+  or  -)  in  Altitude  (km) 

From  Ebq)ected  Values. 

4.  £4  =  0.50;  CH  3  and  CH  4  AT  =  ±5  K  Where:  11 

is  the  Surface  Temperature  (K);  Z^  is 
the  Effective  Cloud  Altitude  (km);  A  £4  is 
the  Difference  (+  or  •)  in  Channel  4 
Emissivity  From  Expected  Values;  and  A  Z 
is  the  Difference  (+  or  -)  in  Altitude  (km) 

From  Expected  Values. 

5.  £4  =  0.90;  CH  3  and  CH  4  AT  =  ±5  K  Where:  12 

T,£.  is  the  Surface  Temperature  (IQ;  Z^  is 

the  Effective  Cloud  Altitude  (km);  A  £4  is 
the  Difference  {+  or  -)  in  Channel  4 
Emissivity  From  Expected  Values;  and  A  Z  is 
the  Difference  (+  or  -)  in  Altitude  (km)  From 
Elxpected  Values. 


1.  INTRODUCTION 

A  cirrus  cloud  analysis  technique  umng  NOAA  Advanced  Very  Hig^  Resolution 
Radiometer  (AVHRR)  multispectral  imagery  has  been  developed  that  determines  the 
bulk  radiative  properties  of  cirrus  clouds.  This  technique,  called  the  AVHRR  Infrared 
Cirrus  model  (AIRC),  computes  optical  depths  and  altitudes  for  ni^ttime  cirrus  cloud 
data  samples.  Cirrus  emissivities  and  altitudes  are  retrieved  using  a  radiative 
transfer  model  that  takes  into  account  both  the  spectrally  dependent  transparent 
nature  of  thin  cirrus  clouds  and  the  attenuative  effects  of  atmospheric  water  vapor. 
The  cirrus  analysis  technique  capitalizes  on  the  differences  between  the  brightness 
temperatures  of  AVHRR  Channels  3,  4,  and  5  (3.7, 10.7,  and  11.8  /tm,  respectively). 
For  scenes  that  contain  thin  cirrus,  the  warmer  underlying  surface  contributes 
significantly  more  energy  to  the  Channel  3  satellite-measured  radiance  than  to  either 
Channel  4  or  5.  These  effects  are  included  in  AIRC  and  allow  for  the  retrieval  of 
cirrus  emissivity  and  hei^t.  The  multispectral  cirrus  analysis  model  demonstrates 
that  more  accurate  thin  cirrus  altitudes  can  be  obtained  than  those  using 
single-channel  longwave  infrared  techniques. 

AIRC  U  restricted  to  nighttime  AVHRR  data  because  incident  solar  radiation 
entering  the  top  of  the  Earth’s  atmosphere  is  not  n^ligible  at  3.7  ^m  and  can,  in 
fact,  be  as  large  as  emitted  terrestrial  radiation,  depending  on  the  temperature 
(Smith  and  Rao,  1972).  It  is  difficult  to  separate  reflected  solar  firom  emissive 
thermal  effects  for  daytime  3.7  /tm  measurements;  this  would  comprise  a  study  in 
itself.  Researchers  at  the  University  of  Utah  are  currently  developing  a  daytime 
cirrus  retrieval  algorithm  (Liou,  personal  communication). 
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A  sensitiyity  analysis  of  the  AIROretrieved  cirrus  parameters  is  mqxulant  to 
users  of  the  technique.  It  provides  insight  into  how  deviations  in  the  correctness  of 
surface  temperature  measurements,  and  also  in  the  accuracy  of  Channel  3  and  4 
brightness  temperature  measurements,  can  lead  to  errors  in  determining  effective 
cirrus  altitude  and  bulk  emissiviiy.  The  methodology  and  results  of  the  sensitivity 
analysis  are  the  focus  of  this  report. 

2.  NOTATION 

B^CT)  Planck  Blackbody  Monochromatic  Radiance, 

Watts  m'*  /Ltm*^  ster'^  [ML'^T®] 

Bj(T)  Planck  Spectrally-Averaged  Blackbody  Monochromatic  Radiance  for 
AVHRR  Channel  j.  Watts  m‘*  /im‘*  ster'^  [ML'^'T*] 

e  Emissivity,  Dimensionless 

Ij  Spectrally-Averaged  Observed  Upwelling  Thermal  Monochromatic 

Radiance  Measured  in  AVHRR  Channel  j.  Watts  m'*  /im**  ster'^  [ML’^T®] 

X  Wavelength,  fim  [L] 

j  AVHRR  Infrared  Channel  Number  (j  =  3,  4,  5) 

r  Reflectivity,  Dimensionless 

I^(X)  Response  Function  at  Wavelength  X  for  AVHRR  Sensor  j,  Dimensionless 
t  Transmissivity,  Dimensionless 
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T 


Temperature,  K  [6] 


X  Atmospheric  Transmittance,  Dimensionless 

z  Heif^t  AGL,  km  [L] 

3.  CIRRUS  PROPERTY  ANALYSIS  BIETHOD 

The  spectrally-averaged  upwelling  thermal  monochromatic  radiance  Ij 
measured  by  a  downward  pointing  radiometer  for  a  field  of  view  filled  completely 
with  a  non-refiective,  thin  cirrus  cloud  is  well  approximated  by 


h  -  *  /  B/7)  ^  <fe 


(la) 


where  j  is  the  AVHRR  infi’ared  (IR)  Channel  number;  e  is  the  cirrus  cloud  bulk 
emissivity  for  AVHRR  Channel  j;  and  ]^(T)  is  defined  by 


where  BJD  is  the  monochromatic  Planck  radiance  at  wavelength  X  for  a  bladdbody 
at  temperature  T;  R](X)  is  the  AVHRR  Channel  j  response  function  within  the 
wavelength  range  s  X  ^  of  sensor  j;  is  the  effective  cirrus  temp^ature;  and 
i^(z)  is  the  spectral  transmittance  for  the  atmospheric  path  between  z  and  the  top  of 
the  atmosphere  at  z  s  «,  defined  by 

kft 

I  KJii.)  dk 

t/d  •  - ! - ,  (ic) 

/  ijPWl  Itfi.)  dX 


where  Xfjiz)  is  the  monochromatic  transmittance  for  the  path  between  level  z  and  the 
top  of  the  atmosphere  at  z  =  »;  is  the  effective  altitude  of  the  cirrus  cloud;  is 

the  ground  emissivity  for  Channel  j;  is  the  ground  skin  temperature;  is  the 
spectral  transmittance  for  the  atmospheric  path  between  level  z  and  the  cirrus  cloud 
base  (defined  similarly  as  tj(z)  above);  and  t^^^  is  the  cirrus  cloud  bulk  transmissivity 
for  AVHRR  Channel  j,  j  =  3,  4,  5. 

The  first  term  on  the  right  side  of  Eq.  (la)  denotes  the  cirrus  cloud  contribution 
to  Ij.  The  second  term  denotes  the  upwelling  emission  by  the  atmosphere  from  the 
top  of  the  cloud  upward.  The  third  term  denotes  that  part  of  the  upwelling  radiance 
emitted  by  the  undertying  surface  at  z=0,  attenuated  the  atmosphere  between  the 
ground  and  the  cloud,  transmitted  through  the  cloud,  and  attenuated  finally  by  the 
atmosphere  above  the  cloud.  The  fourth  term  denotes  the  upwelling  radiance  emitted 
by  the  atmosphere  between  the  groimd  and  the  cloud,  transmitted  through  the  cloud, 
and  attenuated  again  by  the  atmosphere  above  the  cloud. 
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One  important  assumption  that  goes  into  the  fcnmulation  of  Eq.  (1)  is  that  the 
cirrus  cloud  is  non-reflective.  Ice  particle  cloud  radiative  properties  computed  by 
Hunt  (1973)  indicate  that  cirrus  reflectivities  are  veiy  small;  for  an  optically  thick 
cloud  the  Channel  3  (3.7fim)  cirrus  reflectivity  r^  never  exceeds  0.02,  and  Channels 
4  and  5  (10.7  and  ll.S^un,  respectively)  reflectivities  never  reach  0.01.  Thus  for 
cirrus  clouds  the  energy  conservation  equation  is  weU  approximated  by 


or 


fjjhl  “  1  “  ,  y  =  3,  4,  5. 


(2) 


Althou^  the  cirrus  reflectivities  are  effectively  zero  for  Channels  3,  4,  and  5, 
the  emissivities  are  not  as  easily  handled.  Equation  (2)  governs  the  radiative 
properties  of  cirrus  clouds  for  each  channel  individually,  but  offers  no  information  on 
how  the  emissivities  for  a  particular  cirrus  cloud  can  vaiy  from  channel  to  channel. 
Table  1  lists  a  set  of  corresponding  emissivities  for  cirrus  clouds  of  var3ring 
thicknesses,  and  composed  of  spherical  ice  particles  with  a  particle  size  mode  radius 
of  16  fim  (after  Hunt,  1973).  Note  how  the  SJfim  Channel  3  cloud  emissivities  are 
never  larger  than  the  10.7/im  (Channel  4  emissivities.  It  is  possible  to  construct  a 
simple  linear  r^ession  equation  between  the  Channel  3  and  Channel  4 
transmissivities  of  the  form 

fai  =  «  In  +  In  JT  ,  (3) 
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Table  1.  Corresponding  Values  of  Emismvities  for  AVHRR  Channels  3  and  4,  Listed 
as  a  Function  of  Cloud  Thickness.  The  Cloud  is  Assumed  to  be  Composed  of 
SphCTical  Ice  Particles,  and  a  Modified  Gamma  Size  Distribution  is  Used  With  a 
Minimum  Droplet  Radius  of  O.l^m,  a  Maximum  of  25/Am,  and  a  Mode  Radius  of 
l^m;  Number  Density  is  100  cm**  (Aftw  Himt,  1973) 


where  m  and  In  K  are  the  regression  slope  and  intercept  values,  respectively,  in 
1(^-1(^  space.  Performing  the  r^ression  analysis  on  the  transmissivity  pairs 
»  l'®wd)  given  in  Table  1  for  m  and  In  K  and  solving  Eq.  (3)  for  yields 
K=1  and  m=0.67,  or 


=  1  “  (  1 


QjKf 


(4a) 


The  correlation  c  :  '^cient  for  Eq.  (4a)  is  0.997.  This  equation  allows  the  Channel  3 
cirrus  emissivity  to  be  determined  once  the  corresponding  Channel  4  emissivity  is 
known.  The  emissivities  of  Channels  4  and  5  are  more  similar  since  they  represent 
spectral  bands  that  lie  close  to  each  other  in  the  longwave  IR  window.  After  Inoue 
(1987)  it  is  assumed  that  K=1  and  m=1.08,  so  that 

-d (4b) 

which  reproduces  differences  between  Channels  4  and  5  bri^tness  temperatures  in 
the  range  observed  by  Saunders  and  Kriebel  (1988)  for  thin  cimis. 

Consid^  a  pixel  for  which  the  atmospheric  temperature-moisture  profile  and 
underlying  surface  skin  temperature  T^  are  known.  Then  the  transmittance  profiles 
Tj(z)  and  Tj'(z)  can  be  computed  using  a  multivariate  pol3niomial  regression 
transmittance  model  developed  for  AVHRR  th«mal  channels  by  Weinreb  and  Hill 
(1980).  In  fact,  all  terms  on  the  ri^t  side  of  Eq.  (1)  would  then  be  either  known  or 
computable.  Using  Elq.  (2)  to  rewrite  the  cirrus  transmissivity  in  terms  of  the 
emissivily  and  using  Eq.  (4)  to  express  the  Channels  3  and  5  emissivities  as  a 
function  of  the  Channel  4  emissivity  the  satellite-measured  radiance  Ij  can  be 

computed  as  a  function 
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h  ~  «4^ )  1 


y  »  3.  4.  5 


(5a) 


tising  Eq.  (1).  In  Eq.  (5a),  the  upwelling  radiance  ^  is  explicitly  expressed  as  a 
function  of  the  effective  cloud  altitude  and  the  Channel  4  dmis  bulk  emissivity 
Hdd-  observed  radiances  Ij  can  be  converted  to  observed  bri^tness  temp^ature 
Tj  using  a  relationship  similar  to  Elq.  (lb),  so  that 

^  ®  ^  (  t  ^AjfM  )  »  y  ®  3,  4,  5  .  (5b) 


A  sample  set  of  brightness  tenq>erature  calculations  Tg,  T4,  and  Tg  are  plotted 
in  Table  2  as  a  function  of  and  for  an  atmospheric  sounding  taken  at 
Madison,  Wisconsin  on  28  October  1986  at  approximately  0930  UTC.  The  surface 
skin  temperature  T^  was  chosen  to  be  278.4  K  (this  is  1.5  K  lower  than  the  ambient 
air  temperature  of  the  lowest  atmospheric  layer,  a  condition  common  at  ni^t),  and 
the  ground  emissivities  e^,  64^  and  were  chosen  to  be  0.95,  1,  and  1, 
respectively.  Note  in  general  that  for  a  fixed  bri^tness  temperature  measurement 
Tj,  cirrus  altitude  decreases  as  emissivily  increases. 

It  is  possible  to  use  real-time  computations  similar  to  those  summarized  in 
Table  2  in  coqjunction  with  simidtaneously  measured  satellite  brightness 
temperatures  for  a  particular  thin  cirrus  field  of  view  in  order  to  infer  the  cirrus 
effective  cloud  hei^t,  emissivities,  and  transmissivities.  This  is  accomplished  by 
comparing  the  measured  AVHRR  bri^tness  temperatures  to  the  tabulated  ones. 
Consid^  the  Table  2  brightness  temperatures  to  be  representative  of  a  hypothetical 
cirrus  cloud  scene.  Suppose  in  addition  that  the  scene’s  satellite-measured  bri^tness 
temperatures  for  Channels  3,  4,  and  5  are  267.1,  256.4,  and  254.7  K,  respectively. 
Upon  comparison  of  these  bri^tness  temperatures  vdth  the  corresponding 
temperatures  listed  in  Table  2,  the  effective  cloud  hei^t/emissivity  pair 
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simultaneously  consistent  with  the  bri^tness  tenq>erature  measurements  is  =  8 
km  and  =  0.6.  From  Eq.  (4a),  the  corresponding  Channel  3  cirrus  emissivity 
is  0.46,  and  from  Eq.  (4b),  0.63.  Using  Eq.  (2),  the  dmis  transmissivities 

^4,dd>  ^6,dd  0.54,  0.40,  and  0.37,  corresponding  to  optical  depths  and 

^6,dd>  0.62,  0.92,  and  0.99,  respectively  (5j^  •  -In  ). 

Figure  1  contains  a  contouring  of  the  theoretical  bri^tness  temperatures 
(listed  in  Table  2)  as  a  function  of  emissivity  and  effective  cirrus  altitude.  Note  how 
a  Channel  4  brightness  temperature  T4  of  260  K  constrains  the  cirrus  emissivity  to 
be  greater  than  -0.4,  and  the  cirrus  altitude  to  be  hi^er  than  -5  km.  Figure  2 
contains  a  contouring  of  the  theoretical  bri^tness  temperature  differences  listed  in 
Table  2  between  Channels  3  and  4,  plotted  as  a  function  of  emissivity  and  cirrus 
altitude.  Note  how  the  difference  is  at  a  maximum  for  relatively  hi^  values  of 
emissivity  and  altitude  (for  example,  Tj  -  T4  >  26  K  for  an  emissivity  of  0.9  and  an 
altitude  of  12  km). 

For  a  pair  of  coincident  AVHRR  bri^tness  temperature  measurements  Tg  and 
Tj  (j  =  4  or  5 ),  a  set  of  two  equations  in  two  unknowns  and  m  can  be  generated 
using  Eq.  (5b).  Simultaneous  solution  of  these  two  equations  using  a  simple 
graphical  technique  yields  the  (z^j,  E4,dd)  pair  that  corresponds  to  the  measurements 
Tg  and  Tj.  Then  using  Eq.  (4a)  for  j=4  and  (4b)  for  j=;5  the  Channel  3  and  Channel 
5  cirrus  emissivities  can  be  obtained.  Once  Eg^^d  ^j.dd  determined,  Eq.  (2)  can 
be  applied  to  calculate  the  cirrus  transmissivities  tg^^d  ^.dd  approximation 


hjrU 


can  be  used  to  estimate  the  cirrus  biilk  optical  depth  ^j.dd- 
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Table  2.  Channels  3,  4,  and  5  Bri^hti^ss  Temperature  Measurements  for  Various 
Effective  Cloud  Top  Heists  and  Chaund  4  Emissivities,  Calculated  Using  Eqs.  (1) 
and  (4)  for  the  Madison,  Wisconsin  Sanq>le  Atmosphere. 


Effective 

Cloud 

Channel  4  Cirrus  Emissivily 

Height 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

Theoretically-Expected  Bri^tness  Temperatures  Tg  for  AVHRR  Channel  3  (K)  | 

12 

277.5 

274.5 

270.9 

265.9 

257.8 

212.3  1 

10 

277.5 

274.6 

271.0 

266.1 

258.3 

220.8  1 

8 

277.5 

274.8 

271.5 

267.1 

260.3 

237.1 

6 

277.5 

275.3 

272.8 

269.5 

264.9 

253.5 

4 

277.5 

276.3 

274.9 

273.3 

271.2 

267.1 

2 

277.5 

277.6 

277.8 

278.0 

278.2 

278.6 

0 

277.5 

277.8 

278.2 

278.6 

279.0 

279.8 

1  Theoretically-Expected  Brighti^ess  Temperatures  T^  for  AVHRR  Channel  4  (K)  | 

12 

278.6 

269.2 

258.6 

246.5 

231.8 

212.3  1 

10 

278.6 

270.0 

260.4 

249.5 

236.8 

220.8  1 

8 

278.6 

271.8 

264.4 

256.4 

247.4 

237.1 

6 

278.6 

274.1 

269.4 

264.4 

259.1 

253.5 

4 

278.6 

276.4 

274.2 

271.9 

269.5 

267.1  I 

2 

278.6 

278.6 

278.6 

278.6 

278.6 

278.6  1 

0 

278.6 

278.9 

279.2 

279.4 

279.7 

279.9  1 

1  Theoretically-Expected  Brightness  Temperatures  Tg  for  AVHRR  Channel  5  (IQ  I 

12 

278.5 

267.9 

256.4 

243.5 

228.8 

212.3 

10 

278.5 

268.9 

258.4 

247.0 

234.4 

220.8 

8 

278.5 

271.0 

263.1 

254.7 

245.9 

237.1 

6 

278.5 

273.6 

268.6 

263.5 

258.3 

253.4 

4 

278.5 

276.1 

273.8 

271.4 

269.1 

266.9 

2 

278.5 

278.4 

278.4 

278.3 

278.2 

278.1 

0 

278.5 

278.8 

279.0 

279.3 

279.5 

279.7 
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Cirrus  Emissivi ty 

Figure  1.  AVHRR  lO.lfim  Channel  4  Theoretical  Bri^tness  Temperatures  as  a 
Function  of  the  Channel  4  Emissivily  84^  and  the  EfiTectiye  Cirrus  Altitude 


4.  SENSITIVITY  ANALYSIS  METHOD 

The  objective  of  the  remainder  of  this  report  is  to  test  the  AVHRR  infrared 
cirrus  model  for  retrieved  emissivily  cmd  effective  cloud  hei^t  accuracy  when 
introducing  known  ^ors  into  the  independent  variables  used  to  run  the  model. 
Table  2  shows  expected  brightness  temperatures  for  AVHRR  Channels  3, 4,  and  5  in 
K,  for  various  effective  cloud  heights  and  Channel  4  emissivities.  Clearly,  errors  in 
measuring  the  bri^tness  temperature  will  cause  reduced  accuracy  in  the  resultant 
effective  cloud  hei^t  and  cirrus  bulk  emissivity.  As  can  be  seen  in  Figure  1,  the 
range  of  value  in  observed  brightnee  temperature,  over  the  range  of  effective  cloud 
hei^t  hum  0  to  12  km,  increase  e  Channel  4  cirrus  emieivily  increase  from  0  to 
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1.0.  But  this  is  not  a  quantifiable  indication  of  where  the  greatest  error  may  occur, 
'i^gorous  testing  is  necessary  in  order  to  determine  where  the  AIRC  model  is  most 
sensitive.  The  surface  temperature  has  a  significant  influence  on  the  model  as 
well  as  the  aforementioned  AVHRR  brif^tness  tonperature  measurements  T3,  T4, 
and  T5.  The  sensitivity  analysis  takes  into  account  errors  in  T3  and  T4  and  the 
surface  temperature  T^  and  yields  error  ranges  of  emissivity  and  effective  doud 
altitudes.  The  following  error  values  were  examined: 

AT^  =  0,  -t-S,  -3  (K)  (precise,  3  K  too  hi£^  aiui  3  K  too  low,  respectively); 

ATj  =  T3  Measurement  Error  =  0,  +1,  +2,  +3,  +4,  +5,  -1,  -2,  -3,  -4,  -6  (K);  and 


AT4  =  T4  Measurement  Error  =  0,  +1,  +2,  +3,  +4,  +5,  -1,  -2,  -3,  -4,  -6  (K). 


The  true  effective  cirrus  altitudes  were  limited  to  7,  9,  and  11  km  and  to  control  the 
test  size,  Channel  4  emissivities  were  investigated  at  0.2,  0.5,  and  0.9,  ranging  from 
optically  thin  to  nearly  opaque.  Finally,  the  Fort  McCoy,  28  October  1986,  FIRE 
radiosonde  data  were  used  in  running  the  Weinreb  and  Hill  transmittance  model. 
The  results  of  the  sensitivity  analysis  are  presented  in  the  next  section. 

5.  SENSmvnY  ANALYSIS 

The  summaiy  descriptive  statistics  of  the  overall  sensitivity  analysis  are  listed 
in  Tables  3,  4,  and  5.  A  quick  comparison  shows  that  the  largest  errors  occur  at  low 
values  of  emissivity  (£4  =  0.2)  and  dramatically  decrease  in  magnitude  at  hi^ 
emissivity  values  (£4  =  0.9).  At  a  given  emissivity,  induced  errors  in  the  surface 
tCTiperature  have  little  effect  on  the  magnitude  of  the  resultant  deviations  fi*om 
expected  values  of  cloud  altitude  and  emissivity.  These  results  are  useful,  but  to 
obt^  a  better  quantitative  understanding  of  exactly  how  these  variations  are  taking 
place,  we  plotted  the  data  on  equal-area  charts  of  T3  and  T4  errors  AT3  and  AT4, 
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respectively,  as  a  function  of  hei£^t  and  onissivity  (see  Figures  2-10,  Appendix  A). 
The  following  sections  describe  the  infonnation  contained  in  these  plots. 

Table  3.  64  =  0.20;  AT,  and  AT4  =  ±5  K  Where:  is  the  Surface  Temperature  (10; 

is  the  Effective  Cloud  Altitude  (km);  Ae4  is  the  Resulting  Difference  (-i-  or  •)  in 
Channel  4  Emissivily  From  EIxpected  Values;  and  Az  is  the  Resulting  Difference  (+ 
or  -)  in  Altitude  (km)  From  Expected  Values. 


AT* 

*cW 

MAXAz 

AVGAe4 

AVG  Az  1 

(K) 

(km) 

+ 

- 

+ 

- 

+ 

- 

+ 

-  1 

mm 

5.5 

.29 

.16 

3.6 

4.9 

277.1 

M| 

3.5 

.32 

.15 

2.5 

5.7 

CO 

1 

n 

.75 

.20 

1.5 

.32 

.16 

1.4 

7.0 

ig 

13 

5.5 

.26 

3.7 

4.8 

274.1 

3.5 

.31 

^B 

2.2 

5.7 

.S 

II 

0 

.20 

1.5 

.31 

.16 

1.3 

6.8 

5.5 

mm 

.26 

.16 

3.6 

5.6 

271.1 

5.5 

.31 

.15 

3.4 

5.6 

(AT^=-3) 

m 

1.5 

.30 

.16 

1.4 

6.7 
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Table  4.  0.50,  AT,  and  AT4  s  ±5  K  Wh«ne:  T,«,  is  the  Surface  Ten]|>arature  (K); 

Zm  is  the  Effective  Cloud  Altitude  Gan);  AE4  is  the  Resulting  Difference  (+  (u*  -)  in 
Channel  4  Emissiyity  From  Expected  Values;  and  Az  is  the  Resulting  Difforence  (+ 
or  •)  in  Altitude  (km)  From  Expected  Values. 


1  AT* 

Zdd 

MAXAz 

AVGAe4 

AVG  Az  1 

(K) 

(km) 

+ 

- 

+ 

- 

+ 

- 

+ 

ID 

jj^QI 

.49 

.35 

5.5 

R9 

.28 

.21 

3.8 

D: 

277.1 

.45 

.28 

3.5 

.19 

.16 

2.7 

9 

CO 

+ 

i 

.35 

.25 

1.5 

.13 

.14 

1.5 

Bl 

ig 

.46 

.38 

5.5 

2.5 

.23 

1.5 

274.1 

.46 

.29 

3.5 

3.9 

.16 

2.1 

(AT^=0) 

M 

.37 

.26 

1.5 

4.9 

.14 

.15 

n 

2.4 

.50 

.41 

5.5 

.31 

.24 

3.8 

1.5 

271.1 

.48 

.30 

3.5 

.22 

.17 

2.7 

2.0 

(AT^=.3) 

EBSil 

.39 

.26 

1.5 

9 

.14 

.15 

1.5 

i±J 
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Table  6.  £4  =  0.90;  ATs  and  AT4  s  ±5  K  Where:  is  the  Surface  Temperature  (K); 

z^a  is  the  Effective  Cloud  Altitude  (km);  Ae4  is  the  Resulting  Difference  (+<»:-)  in 
Channel  4  Emissivity  From  Expected  Values;  and  Az  is  the  Resulting  Difference  (+ 
or  -)  in  Altitude  (km)  From  Expected  Values. 


AT* 

^cM 

MAXAz 

AVGAe4 

AVGAz  1 

(K) 

(km) 

B 

- 

+ 

+ 

- 

+ 

. 

2.8 

m 

.03 

.09 

0.4 

0.5 

277.1 

.04 

2.2 

.02 

.06 

0.4 

(AT^*+3) 

.02 

.07 

1.5 

n 

.01 

.04 

1.3 

0.5 

.10 

.24 

3.1 

1.0 

.04 

.10 

274.1 

.05 

.12 

3.5 

1.0 

.02 

.06 

99 

(AT^=0) 

SOI 

.02 

.07 

1.5 

1.1 

.01 

.04 

n 

0.5 

.09 

.27 

3.4 

1.0 

.04 

.12 

271.1 

.06 

.13 

3.5 

1.0 

.02 

.06 

(AT^=.3) 

.03 

.07 

1.5 

1.1 

.01 

.04 

BB 

6.  SENSrnvnY  results  for  optically  thin  cirrus  CLOUDS: 
AVHRR  CHANNEL  4  EMISSIVITY  OF  0^ 

6.1  Surface  Temperature  of  277.1  K 

This  section  describes  the  sensitivity  analysis  results  plotted  in  Figure  2, 
Appendix  A. 

When  the  surface  temperature  value  is  in  error  by  +3  K  there  are  large  errors 
in  retrieved  emissivity  (e)  and  cloud  height  (z)  for  low,  medium,  and  hi^  cloud 
altitudes.  The  largest  emissivity  errors  occur  where  the  T3  and  T4  errors  are  both 
n^ative  at  low  altitudes  (7  km)  but  shift  to  where  T,  errors  are  n^ative  and  T4 
errors  are  positive  at  medium  (9  km)  and  hi^  altitudes  (11  km).  In  general,  n^ative 
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errors  in  T3  result  in  a  positive  e  errcr  and  positive  errors  in  T,  result  in  a  n^ative 
e  error  regardless  of  the  magnitude  of  the  T4  error  AT^.  As  the  positive  T,  error 
increases,  the  e  value  decreases  and  goes  to  0.00  for  AT3  2  -t-S  K.  Calculated 
emissivity  values  range  from  0.00  to  0.99  at  7  km  doud  heii^ts  and  0.00  to  0.95  at 
9  km  and  11  km  doud  hei^^ts. 

Cirrus  height  errors  are  present  for  aU  T3  and  T4  error  values  but  become  more 
n^^tive  as  the  true  doud  altitude  increases.  At  7  km,  errors  are  predominantly 
n^^tive  when  both  ATj  and  AT4  are  n^ative,  and  are  predominantly  positive  when 
AT^is  positive  and  AT3  is  n^ative.  At  9  and  11  km,  errors  are  mostly  n^;ative 
when  T3  errors  are  n^ative  for  both  podtive  and  n^ative  errors  of  T4.  Cloud 
altitude  errors  Az^  are  generally  positive  for  0  s  AT3  s  2  K,  and  for  positive  and 
negative  errors  in  T4;  whereas  there  is  no  cirrus  cloud  detected  for  AT3  2  +3  K,  which 
results  in  large  errors.  Calculated  cirrus  hei^t  values  range  from  0  to  12.5  km 
for  all  doud  altitudes  and  errors  increase  in  magnitude  as  the  cloud  altitude 
increases. 

6^  Surface  Temperature  of  274.1  K 

This  section  describes  the  sensitivity  analysis  results  plotted  in  Figure  3, 
Appendix  A. 

If  there  is  no  error  in  the  surface  temp^atme  (that  is,  AT^  =  0  K),  the 
emissivity  errors  that  occur  are  remarkably  consistent  for  all  doud  altitudes  when 
^ors  are  introduced  in  T3  and  T4.  Errors  in  emissivity  are  positive  for  AT3  <  0  K  and 
they  are  negative  for  AT3  >  0  K.  The  largest  positive  e  errors  occur  where  T3  and  T4 
errors  are  both  n^ative  for  a  cloud  altitude  of  7  km;  th^  occur  for  -3  s  3  K  for 
9  km  cloud  altitude;  and  th^  occur  yfhen  AT3  <  0  K  and  AT4  >  0  K  at  11  km  cloud 
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altitude.  Ilie  calculated  e  values  go  to  0.0  for  Ts  arors  exceeding  -f  3  K  The  range 
of  e  errors  is  firom  apfaroodmately  *0.20  to  -f  0.75  for  all  cloud  altitudes. 

Height  errors  at  s  274. 1 K  (that  is,  at  AT^  s  0  K)  are  very  similar  to  hei^t 
errors  at  277.1  K  (AT^  s  -t-S  K).  They  become  increasingly  negative  as  the  true  cloud 
altitude  increases  fi^nn  0  to  12.5  km.  Cloud  hei^t  errors  are  generally  positive  for 
T,  errors  between  0  and  ■l•2  K,  and  f<^  both  positive  and  negative  errors  in  T4.  The 
largest  negative  Az^  occurs  f<»:  AT3  a  h-S  K  since  no  dirus  would  be  detected  at  all. 

6.3  Surface  Temperature  of  271.1  K 

This  section  describes  the  sensitivity  analysis  results  plotted  in  Figure  4, 
Appendix  A. 

Emissivity  errors  are  generally  similar  at  7, 9,  and  11  km  cloud  altitudes  when 
the  surface  temperature  is  in  error  by  *3  K  »  277.1  K);  however  there  are  some 
distinct  differences.  Emissivity  errors  Ae  are  positive  for  n^ative  T,  errors  at  all 
cloud  altitudes  with  the  exception  of  7  km,  which  fails  to  detect  cirrus  at  T4  error 
values  of  +5  K.  Cirrus  is  also  not  detected  for  all  cloud  altitudes  where  AT,  2  -f  3  K. 
Maximum  positive  Ae  errors  occur  at  7  km  doud  altitudes  when  T,  and  T,  errors  are 
both  n^ative  but  shift  to  where  T,  errors  are  n^ative  and  T,  errors  are  positive  at 
9  and  11  km  cloud  altitude.  Emissivily  values  range  from  0.00  to  0.96  for  each  of  the 
cloud  altitudes. 

At  7  km  true  cirrus  altitude,  height  errors  Az^d  are  largest  when  T,  =  +5  K  and 
for  AT,  2  -f  3  this  is  where  the  AIRC  model  would  fail  to  detect  cirrus.  The  model 
would  also  fail  to  detect  cirrus  at  T,  error  values  2  +3  K  for  9  and  11  km  true 
altitudes,  resulting  in  large  negative  height  errors.  As  the  true  cloud  altitude  z^, 
increases  from  7  to  11  km  the  errors  become  more  n^ative  and  increase  in 
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magnitude  firom  -2.8  to  -6.6  K.  The  range  o£  negative  hei^^t  errcnrs  is  from  0  to  -7.0, 
-9.0,  and  -11.0  km  oorreeponding  to  doud  altitudes  of  7,  9,  and  11  km  respectively. 
At  each  altitude  the  errm  values  are  generally  positive  for  T,  oror  values  between 
0  and  +2  K. 

7.  SENSITIVITY  RESULTS  FOR  OPTICALLT  MODERATE  CIRRUS 
CLOUDS:  AVHRR  CHANNEL  4  EMISSIVnY  OF  05 

7.1  Surface  Temperature  of  277.1  K 

This  section  describes  the  sensitivity  analysis  results  plotted  in  Figure  5, 
i^pendix  A. 

For  AT^  s  -t-S  K,  emissivity  errors  are  similar  for  each  of  the  three  true  doud 
altitudes  =  7, 9,  and  11  km.  Positive  onissivity  errors  A64  increase  as  T^  n^ative 
errors  increase,  and  negative  Ae  errors  increase  as  T,  positive  errors  increase.  The 
positive  e  errors  also  increase  as  T4  errors  increase  from  negative  to  positive  values. 
There  is  little  variation  in  n^;ative  e  errors  as  T4  errors  vary,  espedally  at  =  11 
km  altitude.  However,  there  is  a  small  increase  in  n^;ative  e  errors  as  T4  errors 
decrease  from  positive  to  n^;ative  values.  Positive  e  errors  range  from  0.00  to  0.49 
at  2^4  =  7  km,  0.00  to  0.45  at  9  km,  and  0.00  to  0.35  at  11  km;  while  n^ative  e  errors 
have  a  maximum  of  -0.35  at  7  km,  -0.28  at  9  km,  and  -0.25  at  11  km. 

Hei^t  errors  generally  become  more  positive  as  AT,  increases  positively,  and 
become  more  n^;ative  as  AT,  increases  n^atively  for  7, 9,  and  11  km  cloud  altitudes. 
The  hei^t  errors  also  increase  in  a  positive  manner  as  T4  errors  decrease  for  each 
cloud  altitude.  The  magnitude  of  error  for  cirrus  hei^^t  ranges  from  Az^j  s  -3.0  to 
•1.5.5  km  at  7  km;  from  -4.0  to  •f3.5  km  at  9  km;  and  from  -5.0  and  -(■1.5  km  at  11  km- 
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7^  Surface  Temperature  of  274.1  K 


This  section  describes  the  AIRC  sensitivity  anal3rsis  results  plotted  in  Figure 

6,  Appendix  A. 

For  AT^  :=  0  K,  emissivity  errors  at  274.1  K  surface  temperature  are  similar 
to  the  emissivity  otots  discussed  for  AT,fc  =  -3  K  (Section  7.1).  Positive  e  errors 
increase  in  magnitude  as  T,  n^pative  orors  increase  in  magnitude,  and  n^^tive  e 
errors  increase  as  Ts  positive  errors  increase.  The  positive  e  errors  also  increase  as 
T4  errors  increase  while  there  is  little  diange  in  n^ative  e  errors  with  n^pative 
changes  in  T4  oror.  The  magnitude  of  positive  e  mors  ranges  from  0.00  to  0.46  at 
7  km  and  9  km  true  cirrus  altitudes  and  from  0.00  to  0.37  at  11  km  true  cirrus 
altitude.  The  magnitude  of  the  negative  e  errors  decreases  slightly  from  -0.38  at 
=  7  km  cirrus  altitude  to  -0.26  at  =  11  km  cirrus  altitude. 

Hei^t  errors  become  more  positive  as  positive  AT,  increases,  for  all  cirrus 
altitudes  and  for  274.1  K  surface  temperature.  Hei^t  errors  Az^^  become  more 
n^ative  as  positive  T4  error  differencei  increase  at  7  and  9  km  cirrus  altitude  where 
T,  temperature  differences  are  less  than  +4  K  and  at  =  11  km  altitude  where  T, 
temperature  differences  are  less  than  +3  K.  Hei^t  errors  reach  a  maximum  of  +5.5 
km  for  all  AT,  %  4  K  at  7  km  cirrus  altitude;  Az^d  =  +3.5  km  for  all  AT,  2  4  K  at  9  km 
cirrus  altitude;  and  Az^  =  +1.5  km  for  all  AT,  2  +3  K  at  11  km  cirrus  altitude. 
Minimum  height  error  values  range  from  -2.5  km  at  7  km;  -3.9  km  at  9  km;  and  -4.9 
km  at  11  km  cirrus  effective  altitudes. 

7.3  Surface  Temperature  of  271.1  K 

This  section  describes  the  AIRC  sensitivity  analysis  results  plotted  in  Figure 

7,  Appendix  A. 
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Emisuvity  errors  at  271.1  K  surface  temperature  AT^  s  *3  K  are  similar  at  9 
and  11  km  cirrus  altitude  and  are  also  similar  to  the  emissivity  pattern  found  at 
274.1  and  277.1  K  surface  temperature.  For  these  two  altitudes,  emissivity  errors 
become  more  positive  as  Tg  errors  become  more  negative  and  Tg  errors  become  more 
positive.  Emissivity  errors  become  more  n^^ive  as  Tg  errors  increase  positively. 
However,  at  »  7  km  altitudes,  Ac  becomes  more  positive  as  Tg  errors  become  mme 
negative,  and  emissivity  ernurs  become  more  positive  as  Tg  errors  become  more 
positive  for  all  AT  values  except  for  ATg  s  -5  K  and  ATg  s  +5  K.  At  this  point  the 
emissivity  oror  Ac  decreases  instead  of  increasing.  Also,  at  all  values  of  ATg  =  -t-5  K, 
the  Ac  is  positive.  At  values  of  ATg  <  -f  5  K,  e  errors  range  fiem  >0.30  to  +0.48  at  7 
km  dmis  altitude;  from  -0.41  to  +0.48  at  9  km  cirrus  altitude;  and  from  -0.26  to 
+0.39  at  11  km  cirrus  altitude. 

Height  errors  at  9  and  11  km  cirrus  altitude  are  very  similar.  At  2dd=  9  km 
altitude,  Az  errors  decrease  from  a  maximum  of  +3.5  km  to  a  minimum  of  -3.5  km  as 
Tg  error  values  decrease  and  Tg  error  values  increase.  At  11  km  altitude,  Az  errors 
decrease  in  the  same  wiy  as  at  9  km  altitude  but  from  a  maximum  of  +1.5  km  to  a 
minimum  of  -4.7  km.  The  Az  error  is  constant  at  +3.5  km  for  AT,  %  +4  K  and  9  km 
altitude  and  constant  at  +1.5  km  for  ATg  2  +3  K  at  11  km  altitude,  independent  of  Tg 
error  values.  At  7  km  cirrus  altitude,  Az  m'or  values  range  from  a  maximum  of +5.5 
km  to  a  minimum  of  -2.5  km.  The  Az  error  values  become  more  n^ative  as  Tg  errors 
decrease  and  Tg  errors  increase,  with  the  exception  of  one  point  at  ATg  =  -5  K  and 
ATg  =  +5  K  where  the  altitude  error  changes  from  -1.5  km  to  +1.5  km.  The  Az  errors 
do  not  remain  constant  for  any  values  of  AT,  and  ATg  at  7  km  altitude,  as  they 
sometimes  did  for  higher  true  cirrus  altitudes. 
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&  SENSITIVITY  RESULTS  FOR  OPTICALLT  THICK  CIRRUS  CLOUDS: 


CHANNEL  4  EBHSSIVnT  OF  04> 

&1  SorfBoe  Temperature  of  277.1  K 

This  section  describes  the  AIRC  sensitivity  analysis  results  plotted  in  Figure 

8,  Appendix  A. 

For  AT^  =  -t-S  K,  emissivity  error  values  are  very  similar  for  all  altitude  levels 
when  64  =  0.9.  In  general,  positive  e  errors  occur  for  AT3  <  0  K  and  negative  e  errors 
occur  for  ATj  >  0  K.  At  =  7  km,  positive  Ae  values  range  from  0.00  to  0.09;  while 
n^;ative  error  values  range  from  0.00  to  *0.27.  The  largest  negative  e  error  occurs 
as  T4  values  increase  negatively.  Positive  e  errors  increase  sli^tly  as  T4  values 
become  more  positive.  Emissivity  errors  are  almost  symmetrical  at  =  9  km  and 
range  from  0.06  to  0.13.  At  11  km,  the  e  errara  are  qrmmetrical  and  values  range 
from  >0.07  to  0.03,  indicating  a  very  hi|^  level  of  tolerance  for  measurement  error  for 
optically  thick  dmis  clouds. 

Cirrus  hei^t  errors  Az^j  show  strong  similarity  for  each  altitude  examined. 
Hei^t  errors  become  more  positive  as  T4  values  decrease  for  7,  9,  and  11  km 
altitudes.  The  Az  error  values  range  from  -1.0  to  +3.4  km  at  z^a  -  7  km  altitude;  from 
-1.0  to  +3.5  km  at  9  km  altitude;  and  from  -1.1  to  +1.5  km  at  11  km  altitude.  In 
general  about  40%  of  the  Az  error  values  are  n^ative  and  the  remainder  are  either 
zero  or  positive. 

8.2  Surface  Temperature  of  274.1  K 

This  section  describes  the  AIRC  sensitivity  analysis  results  plotted  in  Figure 

9,  Appendix  A. 
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Emisavity  errors  are  small  feu:  AT.^  »  0  K,  and  vary  in  a  similar  manner  as  at 
AT^  s  -fS  K.  Emissivity  errma  vary  from  positive  to  n^^ve  as  AT,  values  go  from 
negative  to  positive  at  7,  9,  and  11  km  true  altitudes.  Emissivity  error  values  range 
fitim  -0.24  to  -fO.lO  at  7  km;  from  *0.12  to  -t-O.OS  at  9  km;  and  fixun  -0.07  to  +0.02  at 
11  km.  The  values  become  more  positive  at  »  7  km  as  AT,  values  become  more 
positive  and  show  little  variation  at  9  and  11  km  altitudes. 

Cirrus  hei|^t  errors  vary  similar^  for  each  of  the  three  altitudes,  in  mudi  the 
same  way  as  the  variation  described  for  277.1  K  surface  temperatures.  The  Az  errors 
become  more  positive  as  AT,  values  decrease  for  each  altitude.  Values  of  Az^,  range 
from  -1.0  to  +3.1  km  at  7  km  true  cirrus  altitude;  from  -1.0  to  +3.5  km  at  9  km 
altitude;  and  from  -1.0  to  +  1.5  km  at  11  km  altitude. 

&3  Surface  Temperature  of  271.1  K 

This  section  describes  the  AIRC  sensitivity  analysis  r^ults  plotted  in  Figure 
10,  Appendix  A. 

Overall,  the  smallest  errors  in  emissivity  occur  at  271.1  K  surface  temperature 
where  e,  =  0.9.  Values  of  Ae  range  from  -0.20  to  +0.09  at  7  km  true  cirrus  altitude; 
from  -0.10  to  +  0.04  at  9  km  altitude;  and  from  -0.07  to  +0.02  at  11  km  altitude.  At 
=  7  km,  there  is  a  small  positive  increase  in  emissivity  error  values  as  T, 
temperature  errors  become  more  positive,  but  very  little  change  in  e  error  values 
occurs  as  AT,  values  change  at  both  9  and  11  km.  Again,  values  of  e  error  are 
positive  where  AT,  values  are  n^;ative  and  are  negative  where  AT,  values  are 
positive. 

Overall,  the  hei^t  error  values  Az^  are  also  the  smallest  for  271.1  K  surface 
temperature  and  e,  =  0.9.  Values  of  Az^,  range  from  -1.0  to  +2.8  km  at  z^,  =  7  km 
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true  altitude:  -0.9  to  +2.2  km  at  s  9  km  altitude;  and  from  -1.1  to  +1.5  km 

at  ZcM  s  11  km  altitude.  In  general,  Az  error  values  are  n^^ative  where  AT^  values 
are  positive  and  are  positive  where  AT4  values  are  n^ative.  The  variation  in  Az  error 
values  is  very  similar  to  those  discussed  earlier  for  £4  =  0.9  at  277.1  K  and  274.1  K 
surface  temperatures. 

9.  SUMBIARY 

In  general,  retrieved  emissivity  and  altitude  errors  are  greatest  when  the  true 
cirrus  emissivity  £4  is  low,  and  least  when  £4  is  hi^.  Therefore,  the  AIRC  cirrus 
model  is  least  sensitive  to  T3  and  T4  measurement  errors  when  the  cirrus  is  optically 
thick.  Errors  in  both  effective  emissivity  and  altitude  (Ae  and  Az^^  respectively)  are 
greatest  at  £4  =  0.2  (optically  thin  cirrus).  In  all  cases  the  errors  decrease  as  the  true 
altitude  z.^^  increases  from  7  km  to  11  km.  The  errors  also  tend  to  decrease  as  the 
surface  temperature  values  increase  from  -3  K  to  +3  K  of  the  true  T^  value.  Finally, 
less  variation  in  error  occurs  for  hi^er  than  for  lower  £4  values  and  in  general  the 
largest  error  values  occur  at  £4  s  0.2. 

These  error  tendencies  are,  in  part,  due  to  the  constraints  placed  on  cirrus 
attributes  by  AIRC.  The  decrease  in  altitude  error  from  7  to  11  km  is  affected  by  tiie 
limitation  of  a  12.5  km  height  for  cirrus  clouds  imposed  by  the  model.  Also,  it  is 
logical  for  the  largest  errors  to  occur  for  the  optically  thinnest  cirrus  (£4  =  0.2),  since 
the  measured  AVHRR  bri^tness  temperatures  are  minimally  affected  by  the  cloud 
itself. 


Somewhat  surprisingly,  the  accuracy  of  the  surface  skin  temperature  T^  turns 
out  to  be  one  of  the  less  sensitive  aspects  of  the  AIRC  cirrus  model. 

Surface  emissivities  can  be  veiy  difficult  to  specify,  especially  for  the  3.7 
fim  Channel  3.  Aside  from  a  dominant  wavelength  dependence,  varies  on  local 
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and  even  subpixel  scales  due  to  efifects  such  as  variable  soil  moisture  content, 
urban/rural  mixed  fields  of  view,  v^tation  cover,  and  the  like.  Cirrus  cloud 
emissivities  and  transmissivities  have  complex  dependencies  on  wavelength,  particle 
size  distribution,  and  ice  particle  riiapes.  Errors  in  the  vertical  temperature  profile 
or  the  atmospheric  transmittances  will  manifest  thmnselves  largely  in  the  accura<7 
of  the  retrieved  cloud  heights;  the  effect  on  dear-column  radiance  conq)utations  is 
small  since  each  of  the  three  AVHRR  IR  channels  lie  within  reasonably  clean 
atmospheric  windows. 

Sources  of  error  in  measured  bri^tness  temperature  can  stem  fi*om  the  use 
of  AVHRR  local-area  coverage  (LAC)  1  km  data.  LAC  footprints  fi^m  one  AVHRR 
channel  to  the  next  are  usually  not  predsely  aligned,  meaning  that  the  "coinddent" 
brightness  temperature  measurements  T^,  T4,  and  T5  mi^t  not  be  from  precisely  the 
same  FOV.  One  way  to  alleviate  this  problem  is  to  average  a  set  of  3  x  3  or  4  x  4 
LAC  pixels  together  before  inserting  them  into  the  AIRC  cirrus  algorithm.  Care 
would  have  to  be  exerdsed,  however,  to  ensure  that  each  of  the  pixels  being  averaged 
is  completely  filled  with  cirrus.  (This  is  an  important  assumption  of  the  AIRC 
model.)  Finally,  Channel  3  is  known  to  be  somewhat  noi^;  again,  this  problem  may 
be  alleviated  through  the  application  of  spatial  averaging  techniques. 

Spedal  thanks  is  given  to  Dr.  Michael  K.  Griffin  (PL/GPAS)  who  provided  the 
program  used  in  generating  the  plots  of  ^pendix  A.  The  most  ^ective  use  of  this 
report  is  for  AIRC  dmis  model  usera  to  refer  to  these  plots  for  a  perception  of 
possible  errors  inherent  in  the  model  removal  results. 
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APPENDIX  A 

PLOTS  OF  AIRC  SENSITTYITY  ANALYSIS  RESULTS  AS  A 
FUNCTION  OF  T,  AND  T^  MEASUREMENT  ERRORS 
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Figure  3.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  AT,  and  AT4.  Error  Plots  are  for  T,,.  ^  274.1 
K  (AT,fc  =  0  K)  and  True  Cirrus  Emissivity  =  0.2. 


EMISSIVITY  ERRORS  EMISSIVITY  ERRORS  EMISSIVITY  ERRORS 

T-sfc  =  271.1  K  e  =  0.2  z  ■  7.0  km  T-sfc  -  271.1  K  e  -  0.2  *  =  9.0  km  T-sfc  »  271.1  K  e  *  0.2  *  »  11. 
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Figure  4.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  AT,  and  AT4.  Error  Plots  are  for  =271.1 
K  (AT,fc=  -3  K)  and  True  Cirrus  Emissivity  =  0.2. 


EMISSIVITY  ERRORS  EMISSIVITY  ERRORS  EMISSIVITY  ERRORS 

T-sfc  =  277.1  K  e  =  0.5  z  =  7.0  km  T-sfc  -  277.1  K  e  »  0.5  z  .  9.0  km  T-sfc  =  277.1  K  e  >  0.5  z  .  11.0  km 
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Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  ATg  and  AT4.  Error  Plots  are  for  T.^,  =  277.1 
K  =  +3  K)  and  True  Cirrus  Emissivity  =  0.6. 


Figure  6.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Bri^tness  Temperature  Measurement  Errors  AT,  and  AT^.  Error  Plots  are  for  T.^  .  274.1 
K  0  K)  and  True  Cirrus  Emissivity  =  0.6. 
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Figure  7.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  AT,  and  AT4.  Error  Plots  are  for  T,|^  s  271.1 
K  (AT,fc=  -3  K)  and  True  Cirrus  Emissivity  64  «  0.6. 
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EMISSIVITY  ERRORS  EMISSIVITY  ERRORS  EMISSIVITY  ERRCMIS 

T-sfc  =  274.1  K  e  =  0.9  z  =  7.0  km  T-sfc  =  274.1  K  e  =  0.9  z  =  9.0  km  T-sfc  •  274.1  K  e  .  0.9  z  -  11.0  km 
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Figure  9.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Hei^t  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  AT,  and  AT4.  Error  Plots  are  for  T.^,  *  274.1 
K  (AT,fc=  0  K)  and  True  Cirrus  Emissivity  64  =:  0.9. 
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Figure  10.  Plots  of  AIRC  -  Retrieved  Cirrus  Emissivity  Errors  (dimensionless)  and  Cloud  Height  Errors  (km)  as  a 
Function  of  Brightness  Temperature  Measurement  Errors  ATj  and  AT4.  Error  Plots  are  for  T.^  =  271.1 
K  (AT,fc=  -3  K)  and  True  Cirrus  Emissivity  =  0.9. 


